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Abstract. Cosmo logical inhomogeneities gravitationally deflect radiation propagating 
from distant sources, transforming the spatial and angular correlation functions of intrin- 
sic source properties. For a gaussian distribution of deflections (e.g. from a primordial 
gaussian density perturbation spectrum or from the central limit theorem) we calculate 
the probability distributions for geodesic deviations. If the intrinsic variable is also gaus- 
sian, e.g. the large scale velocity flow field or cosmic microwave background temperature 
anisotropies, then distributions and correlation functions of the observed image sky prop- 
erties can be obtained. Specialising to CMB temperature fluctuations we rederive simply 
the influence of independent gravitational lensing on the anisotropy angular correlation 
function and calculate the new effect of lensing correlated with the anisotropies, e.g. aris- 
ing from the same primordial gravitational perturbation field. Characteristic magnitudes 
and scales are given in terms of the density power spectrum. The correlated deflection- 
temperature effect is shown to be negligible. 
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1. Introduction 

The origin and development of structure in the universe needs to be probed by means of 
observations of distant sources. Not only the direct properties of those sources, however, 
but also the characteristics induced in the radiation propagating through the gravitational 
inhomogeneities of the structure provide clues to the nature and formation of structure. 
Within linear theory, i.e. when the deviations from homogeneity are still perturbations, 
there are simple and intimate relations between the important cosmological descriptors of 
the matter, such as the large scale velocity flow field or the energy density fluctuations, and 
of radiation propagation characteristics such as temperature anisotropies (redshifts) and 
gravitational lensing deflections (momentum changes). These interrelationships constrain 
theories of structure formation. 

This paper concentrates on the statistical effects of the inhomogeneities on obser- 
vations, developing an analytic approach that requires only that the physical processes 
behind the observable and lensing quantities be gaussian in nature. This avoids sensitiv- 
ity of results to a specific geometry or density perturbation model and is both physically 
relevant and surprisingly powerful. Although here we are primarily concerned with tem- 
perature anisotropies in the cosmic microwave background radiation (CMB), the method 
is easily generalised to investigations of the velocity flow field or density perturbations. 

Before primordial CMB anisotropies were first detected by the Cosmic Background 
Explorer satellite (COBE; Smoot et al. 1992) there was much discussion in the literature 
concerning the effects of inhomogeneities and their gravitational lensing on the anisotropies 
(Blandford & Narayan 1992 provide a review with references). More recently Fukushige, 
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Makino, and Ebisuzaki (1994) investigated this by calculating scattering with an N-body 
code, and Tomita (1996) considered the effects of superhorizon scale inhomogeneities. For 
the statistical distribution of the lensing deflections, in particular the correlation of neigh- 
boring lines of sight, crucial to the key property of geodesic deviation, various assumption 
have been made, ranging from a diffusion approximation to nearest neighbor correlations. 

No previous work, however, dealt with the case where the deflections are not super- 
imposed independently on the temperature fluctuations but are physically correlated with 
them. For example, consider a hot spot in the CMB. This corresponds to a region where 
the primordial gravitational potential deviates strongly from the mean but that potential 
also determines the strength of the lensing caused by the primordial density field from 
that region. Thus one can imagine that extremes of the temperature field are preferen- 
tially strongly lensed relative to the milder deviations, leading to a significant distortion 
of the intrinsic temperature correlation function. 

This paper investigates the properties of cosmologically important gaussian fields in- 
cluding the possibility of such a cross correlation. The joint probability formalism allows 
for concise derivations of familiar results, e.g. beam smearing, from gravitational lensing 
while showing new effects from the cross correlation. By estimating the magnitude of these 
on anisotropy observations we can probe large scale structure in the universe by placing 
constraints on the density power spectrum The conditional probability formalism 
also points out ways to test the underlying gaussian nature of the primordial cosmological 
perturbations. 

In §2 we review the mathematical basis and derive properties of joint and conditional 
probability distributions for the variables entering the problem - the relative deflections 
of the null geodesies and the intrinsic source sky characteristics. These are combined 
in §3 to form expressions for the observable correlation functions, especially under the 
condition of coherence between the source and propagation conditions. Section 4 relates 
the mathematical results to the underlying physics in terms of the density fluctuation power 
spectrum and §5 presents quantitatively the effects on CMB anisotropy measurements. 

2. Probability distributions 

— * 

In an imperfectly homogeneous universe radiation observed at some angular position 
on the sky may have been deflected during propagation to us from a source whose true 
position is ip. By true position we mean that position at which the source would appear if 
the inhomogeneity (gravitational lens) were smoothed out. A general line of sight has 



for a sky projected deflection angle 9, with 9 varying along different directions. The central 
question for observations in a universe with inhomogeneities is how to relate an intrinsic 
function of the source position to the observed function of image position. 

In order to investigate this in as model independent a way as possible, we assume that 
the gravitational lensing deflection and the physical source variables of interest individually 
have d dimensional gaussian distributions with zero mean (no preferred direction). Then 



$= * + 9 



(1) 



p(v) = [n(2a 2 /d)]- d / 2 e- v2 /^ 2 / d \ 



(2) 
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where the variance a 2 = (v 2 ). One justification for adopting a gaussian is that if the 
variable, e.g. projected deflection 9, is compounded out of many elements, e.g. isolated 
spatial deflections at different distances from the observer (many inhomogeneities along 
the line of sight), then the central limit theorem leads to a gaussian distribution. Another 
rationale enters if the underlying physical process has a gaussian nature, generated by 
quantum fluctuations for example. One case would be gaussian density perturbations which 
create gaussian distributions in lensing deflection, velocity fluctuations, and temperature 
anisotropies in the linear regime. This relation is discussed further in §4. 

Since absolute shifts in the source position are undetectable we must compare at least 
two light rays, which requires knowing the joint probability of deflections along lines of 
sight 1,2,.... More generally one can consider the variables Vi to denote different physical 
processes, such as density or temperature fields. Using the characteristic function approach 
(Cramer 1957) the joint probability of occurence is found to be 

p(i?i, *&,...,*„) = (27r/rf)-^/ 2 M- d / 2 exp{-A^M^-^}, (3) 

where M and My are the determinant and cofactors, respectively, of the covariance matrix 
composed of elements rriij = (vi ■ Vj). Angle brackets denote integration over the joint 
probability distribution. 

In the case of two sky variables (n = d = 2), for example comparing deflections of two 
rays, (3) simplifies to 

p0 u &) = [Tr 2 ^^^! - 6 2 ) j-i e -[(^/ CT ?)-26(e 1 -^ 2 /-i- 2 )+(^/-g)]/(i-fo 2 ) 7 (4) 

where the correlation coefficient b = (9\ • #2)/ciC2- This distribution is normalized, and 
integration over one variable recovers the gaussian distribution for the other variable. As 
6^0 (incoherence or unrelated variables) the distribution separates into two independent 
gaussians while as b — > 1 (total coherence) it approaches a delta function 5(9% — 9\). 

One can imagine physical situations in which the variances <ti , 02 of the two variables 
6*i, 6*2 are not identical even though they correspond to the same physical process, e.g. 
gravitational deflection. For example if the universe were not homogeneous then lines of 
sight through different patches of the universe might have different statistical properties. 
Even with homogeneity, consider sources at different redshifts. Because of the differing 
path lengths the cumulative gravitational deflections and hence variances a\ , 01 would be 
unequal, as is evident from the explicit physical expressions (34) and (36). If we concentrate 
on the source being the CMB, however, the variances of deflections along any lines of sight 
are equal: (9f) = (9 2 ) = cr 2 . 

For other physical situations and gaussian processes, such as the large scale velocity 
field, evaluation out to different distances can be of interest, e.g. windowing different 
volumes of the density field or slicing redshift surveys. Then, for example, one can find 
probability distribution functions such as for the ratio of the variable value in sample 1 to 
that in sample 2, k = (0i/cri)/(0 2 /c2): 

p(k) = 2(1 - b 2 ) k~\k + k- 1 ) [(k + k' 1 ) 2 - 46 2 ]" 3 / 2 . (5) 
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As expected, the differential probability p(k) dk is invariant under transforming k to its 
reciprocal. One can also calculate angular distributions, e.g. p(fj), where fx = 9\ • 62/6162- 
But this is of more use when dealing with three dimensional vector variables such as vr, 
the velocity flow field out to some survey depth R. 

Of particular physical interest is the difference between a quantity evaluated along 
two different lines of sight, e.g. the relative deflection or deviation D = 62 — 6\. Similarly 
one can investigate the coherence of the variables, e.g. (3 = 6\ • 62- Both are important to 
the physics. From (1) the observed separation of two rays $ = <p2 — 4>i is 

$ = tf+D, (6) 

where ip = ^2 — ^1 is the intrinsic (true) separation and the separations can represent 
either the angular distance between two sources or the distance within a source (e.g. source 
size). Thus D is involved in the transformation of field separations and the quantities that 
depend on them. The coherence variable (3 measures how independent the deflections are. 
Physically, we expect that neighboring rays have strongly correlated deflections while lines 
of sight far apart have independent deflections. 

To find the statistical distribution of the deviation and coherence one changes variables 
in (4) to, say, D, /?, and h = 6\ — 9f, employing the Jacobian of the transformation: 

d 4 P(D, p, h) = J- 1 (^M) P(6i, 62) 

\ 61,62 J (7) 

= (nae^Q-'e-^^e-^-^^iD 4 + A(5D 2 - h 2 )- 1 ' 2 d 2 Dd(3dh, 

with Q = oJ(l ~ b 2 ). From this new joint distribution one derives the individual distribu- 
tions by integrating over the other variables: 

p(D) = [27ra|(l - 6)]" 1 e-^ 2 / 2CT ^ 1 - 6 ), (8) 

p(h) = [2al^l^Y 1 e-^l"i-^ rS ' . (10) 

Equations (8-10) illustrate that the distribution of a linear combination of gaussian vari- 
ables is itself gaussian (Cramer 1957); note only D is such a linear function. The magnitude, 
or modulus, of a gaussian distributed vector quantity is Rayleigh distributed, 

p(D) = 2(D 2 )- 1 De- D2 ^ D2 l (11) 

We can find the moments of our chosen variables: 

(D) = = (D 2n+1 ) ; (D 2 ) = 2aj(l — b) ; (D 2n ) = n\ (D 2 ) n , (12) 
((3)=ba 2 9 ; ((3 n ) = 2~ in+1) nl* 2 e n [(1 + 6) n+1 + (-) n (l - 6) n+1 ] , (13) 
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and verify that they have the expected behavior in the limits 6 = and 6=1. The above 
use of relative and coherence variables, Jacobians, and moments is equally applicable to 
variables besides the ray deflection. 

Note that since by (6) the observed separation on the sky depends only on the intrinsic 

— * 

separation and D, we can write correlations of position dependent quantities using only 
p(D) and not the full p(6>i,6> 2 ). For example for any function / only of the relative field 
positions, e.g. separation (which is all that can enter in a globally homogeneous universe), 



J d% J d%f(e 2 -e 1 ) P (e 1 ,e 2 ) = j "d 2 Xj d 2 D f(3) P (D,A) 

= J d 2 Df(D)p(D), 



(14) 



where, say, A = (6\ + 2 )/2. 

We can immediately see that for the transformation of an intrinsic sky function Co of 
only the source separation ip into an observed function C of the image separation <E>, 

C($) = J d 2 D J d 2 ifC (ij)p(D)5($>-iJ-D) 

= ! d 2 $C Q {il>)p{D = $-$) 

r , , (15) 

= (7rS 2 )- 1 / d 2 ^C {^)e-^-^' s2 

POO 

= 2S~ 2 / #^ (f)e-^ +$2 )/ s2 /o(2#/5 2 ), 
Jo 

where Iq is a modified Bessel function and S 2 = (D 2 ). 

This recreates equation (24) of Wilson & Silk (1981) for beam smearing, or scattering 
due to the observer and hence independent of source and path properties. In this case S 
is just a constant corresponding to the dispersion of the telescope beam response. The 
formalism epitomized by (15) also provides a compact method for deriving equation (19) 
of Linder (1990a) for the independent lensing of the microwave background temperature 
anisotropy correlation function. Here the path properties do enter but are uncorrelated 
with the intrinsic source characteristics observed. Then S will be a function of path length 
and field separation due to the different physical conditions along different lines of sight. 

The peak value of the correlation function, at $ = 0, will be diluted by a factor 
1 + 5' 2 /(20 2 ), where G is the characteristic scale, or coherence angle, of the observed 
quantity (cf. equation 18). Any oscillations in the angular behavior Co(ip) will be damped 
by the coarse graining of the integral. Observers seek to have beam sizes much smaller 
than the anisotropy scale of interest, S G, in order to avoid this smearing. In the 
independent lensing case the zero lag correlation is unaffected since at $ = the lines of 
sight are identical, forcing S = (6 = 1). Smearing does enter at larger angles though 
(cf. Linder 1990a, Figures 2 and 3), its importance going as S 2 /Q 2 ~ (ae/Q) 2 (l — 6). 

However, (15) only holds if the lensing process is independent of the observed process, 
i.e. detected characteristic, so that the separation into Co(ip) p(D) is possible. Section 3 
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provides more detailed discussion of this point and treats the correlated case, correcting 
(15) to (24), (25). Meanwhile let us continue examining the necessary variables. 

The other gaussian variable to enter our correlated microwave background calculation 
is the fractional temperature anisotropy U = AT(<f>i)/T. This is a scalar variable (d = 1) 
on the sky. Hence its probability function is 

p(t u t 2 ) = (2 7 r ( 7 t 2 v / r :: ^)- 1 e- (t? - 2Btlt2+ ^ )/2CTt2(1 - B2) , (16) 

where the variance a 2 = (t 2 ) and the correlation parameter B = (tit 2 )/a 2 . In a more 
explicit notation, 

/ AT ^ AT ^ \ 
^ 2 =(— (00— (0 2 )) =<?(*), (17) 

\ 1 1 I *=|£i-<? 2 | 

the familiar angular anisotropy correlation function. In this notation, 

a 2 =(t 2 ) = C(0) 

B = (ht 2 )/a 2 = C($)/C(0) (18) 

e = [-c"(o)/c(o)]- 1 / 2 = (-d 2 B(o)/d$ 2 )- 1/2 , 

relating the correlation function C and the coherence angle or curvature © to the prob- 
ability function characteristics a t and B. Because the temperature is a random process 
evaluated at the same last scattering surface regardless of line of sight, the variance a 2 is 
not dependent on direction (cf. equation 31). 

While B or C($) gives the temperature anisotropy correlations averaged over the 
sky, i.e. the rms value, one can also calculate the entire probability distribution of the 
correlations from (16). For the coherence variable c = tit 2 this is 

p(c) = (Tr^v^TBaj-^oIlcl/af (1 - B 2 )] g^d-s 3 ), (19) 

where Kq is a modified Bessel function. This is important for computing the cosmic 
variance or intrinsic error due to the fact that we observe only a single universe, i.e. one 
realization of the probability distribution. While the mean (c) = o 2 t B = C, the variance is 

a 2 c = (c 2 )-(c) 2 = at(l + B 2 ). (20) 

Note that the distribution of c is not gaussian. 

Observations of the CMB sky temperature can be related to C according to the 
geometry of the fields, i.e. the number of telescope beams and the angles between them. 
This can be made explicit by defining the differential variable r = t 2 — t\ with gaussian 
distribution 

p( r ) = [47ra 2 (l - S )]-V2 e -r 2 /^ 2 (i-B)_ (21) 
A two beam experiment with throw $ measures the variance 

(r 2 ) = 2a 2 (l - B») = 2 [C(0) - C(d>)], (22) 
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while a three beam one symmetric about the center field measures the variance 



(R 2 ) = ((r a -r b ) 2 ) = 2(r 2 )(l-b ab ) 

= 2a 2 (3 - + S 2$ ) = 6C(0) - 8C($) + 2C(2$), 

where indices a, 6 represent the two pairs of fields. Note the similarity of form between 
(22) and (23): a\ = (t 2 ) -> (r 2 >, Ba t 2 = fafe) - (r a r 6 ) = 6 a6 (r 2 ). 

Rather than obtaining only the variance one could, if desired, compute other quantities 
such as the conditional probability of obtaining a certain value in one set of fields given that 
a fixed value was observed in other fields. This could be useful in that discrepancies between 
observations and probability theory could be a possible probe of nongaussianness in the 
temperature fluctuation distribution since only gaussian processes are fully determined by 
the variance. 

In summary, for the physical situation of gravitational inhomogeneity deflection of mi- 
crowave background radiation we are interested in the mathematical correlations between 
two two dimensional gaussian variables - one (the sky deflection 9) a vector and one (the 
temperature fluctuations t) a scalar. We will see that it is in fact convenient to decompose 
the vector deflection into arbitrarily oriented but orthogonal components, thus giving a 
total of three variables per field, all magnitudes. 

3. Cross correlations 

The expression for the observed CMB temperature anisotropy correlation function will 
involve both the intrinsic temperature correlations and the conditional probability distri- 
bution of deflections given that the intrinsic fields have certain temperature fluctuations, 
i.e. values of the gravitational potential field. 

Mathematically, the intrinsic correlation function is defined by the integral of the field 
temperatures weighted by the joint probability of obtaining those particular values in each 
field, 

C (VO = J dti J dtzpih^htz, (24) 

— * 

where ip is the field separation. To include the gravitational deflections one simply maps the 
intrinsic sky coordinates onto the observed ones because quasistatic gravitational lensing 
induces no redshift, i.e. the photon energy and hence temperature is unaffected. The 
geometry of this mapping is illustrated in Figure 1. 

Now to calculate the observed correlation function one employs Bayes' theorem and 
simply adds up all the ways in which the fixed observed field separation $ can be gen- 
erated from the various intrinsic field separations ip. That is, one sums over all possible 
quadrilaterals with one fixed side, weighted by the probability for achieving that shape, 
i.e. obtaining deflections of the necessary magnitudes and directions, given the physical 
conditions represented by t\,t 2 - 

So the observed correlation function becomes 

C($) = J d 2 D J d 2 i) J dh J dt 2 p(t u t 2 ) t 1 t 2 p(D\t u t 2 ) 5 (2) ($ -f-D), (25) 
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where D = 9 2 — 9\ is the relative deflection mapping ip — > <&, the delta function enforces 
closure of the quadrilateral, and p(D\ti,t2) is the conditional probability. This is a direct 
translation of the above words into an equation. For example the phrase "given the physical 
conditions" indicates that the probability function of all the variables does not range over 
the whole parameter space but is constrained by some prior knowledge or condition; hence 
one uses the conditional probability distribution to reflect this. Now conditional probability 
is calculated by the full probability divided by the probability to achieve the set condition. 
Hence we can convert the product of the temperature probability function p(t±,t2) and 
the conditional deflection probability function given those temperature values p(D\ti,t 2 ) 
to a joint probability over all the variables p(D, ti, t 2 ). 

However, the expression is still not that useful for comparison with observations as 
it requires an integration over the ensemble of probability realizations, yet we have only 
a single universe to work with. The solution is to employ the ergodic hypothesis which 
states that we can substitute an integration over the whole sky for an ensemble average. 
Thus we can replace integration over all the values which the random process could take 
on at a location by integration over all evaluation locations. 

The final form of the observed correlation function is then 

C($) = / dV / dt i / dt 2 t 1 t 2 p(D = ^ - tf,t 1 ,t 2 ) 

J J J (26) 

= / dV / d 2 if 1 t(if 1 )t(if 1 + if)p(D = $-if,t 1 ,t2), 

enforcing ip 2 = — ipi in the second line. Only in the case where the deflection process 
is independent of the temperature fluctuations does the joint probability separate into 
a product of p(D) p(ti, £2), necessary to obtain the explicit appearance of the intrinsic 
correlation function Co(V') in (15). 

To obtain the joint probability function between the gaussian variables of the ray 
deflections and the sky temperatures, it is convenient to decompose D into orthogonal 
components. This both puts all variables on an equal dimensional footing and makes the 
correlation (D x D y ) = 0. By using the relative variable D = 62 — #1 we no longer explicitly 
need the deflection correlation parameter f3 = 6\ -6 2- In §4 we find that only the component 
of D along the field separation has a nonvanishing correlation with the field temperatures, 
so we choose our decomposition parallel and perpendicular to this axis. 

Now the joint probability function for the variables D\\, D±, £1, t 2 follows from (3) via 
the covariance matrix, 



M(D h D ±1 t u t 2 ) = 



/ s 2 -bsa bsa \ 

s 2 

-bsa a 2 Ba 2 

V bsa Ba 2 a 2 ) 



(27) 



One can read off from the matrix that s 2 is the variance of each component of the deflection 
.D, a 2 is now the variance of the temperature fluctuations, B is the temperature corre- 
lation parameter, and b is now the cross correlation parameter between the gravitational 
deflection and the temperature fluctuation. 
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The joint probability distribution is then 



p(D h D ± , h, t 2 ) = (2ivsa)- 2 (l + B)~ 1/2 (l -B- 2b 2 )~ 1/2 x 

exp{-£/[2sV(l + B)(1-B- 2b 2 )}}, (28) 
E = s 2 (l - b 2 ) (t\ + 4) - 2s 2 (B + b 2 ) tit 2 - 2bsa(l + B) D\\ (t 2 - h) 
+ a 2 (l - B 2 ) D\ + a 2 (l + B)(l - B - 2b 2 ) D 2 ± . 

Carrying out the integrations specified in (26) yields the observed CMB temperature 
anisotropy correlation function 

„2 



C($) = (2ns 2 )- 1 J d 2 ^C {^)e-^-^ 2 / 2s 

+ (2tts 2 )- 1 J d 2 $b 2 C (0) [l-{($-^)-^} 2 /« 



e _|$_^|2 /2s 2 



— Cind ~\~ Ccor- (29) 

The first term is proportional to B = Co(^)/c 2 and is the contribution to the observed 
temperature correlations when the gravitational lensing is independent of the temperature 
anisotropies, i.e. there is no cross correlation. This is precisely the old result derived 
in Linder (1990a) and mentioned in (15). The second term, involving 6, represents the 
new effect of including correlations between the gravitational deflection field and the tem- 
perature fluctuation field, for example due to both originating from the same primordial 
gravitational potential perturbation. 

The two crucial ingredients are thus the correlations B [or Co(V0] an ^ The former is 
a scalar function of field separation, i.e. independent of orientation, and from its definition 
one expects the latter to vary with angle as 3> • ip; both are verified in the next section. 
Then the angular part of the integrals can be carried out without knowing their explicit 
form: 

poo 

C($) =s- 2 e-* 2 / 2s2 / di;i;[CoWlo(^/s 2 ) + Co(0)P 2 ^^)Y}e-^/ 2s2 

Y = [($ _1 + 2$s -2 ) V + (2$ + 4s 2 $ _1 ) if;' 1 + 6s 4 $ V" 3 ] h{^/s 2 ) ^> 
- [s~V + (1 + $V 2 ) + 3s 2 V>" 2 ] Io($ip/s 2 ), 

where I v are modified Bessel functions and (3 = &/(<& • ip). 

The prescriptions for the variable variances and correlation parameters are given by the 
physics underlying the random processes and can be written in terms of the gravitational 
potential perturbations, or more conventionally the primordial density fluctuation power 
spectrum, as in the next section. 

4. Density fluctuation spectrum 

The mathematical basis of the statistical moment and correlation analysis relied on the 
gaussian nature of the observational characteristics. This implied that all the physics 
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resided in the variance and cross correlation of the characteristic variables. Now we inves- 
tigate those physical processes to obtain expressions for these quantities in terms of the 
underlying mechanism - fluctuations in the cosmological density field - expressed in terms 
of the density power spectrum. 

First consider the source property of the cosmic microwave radiation background tem- 
perature anisotropies. The temperature fluctuations, or shifts in the photon energies, are 
induced by variations in the gravitational potential at the last scattering surface, t = 0/3, 
(variations along the photon path contribute negligibly due to their quasistatic nature). 
In turn, the potential can be related to the density perturbations via the Poisson equa- 
tion, reading (fi k = (3/2)Hq k~ 2 5 k in Fourier space. Thus the variance of the temperature 
fluctuations can be written in terms of the density power spectrum P k = \5 k \ 2 as 

poo 

a 2 = (t 2 ) = C (0) = (Stt 2 )" 1 ^ / dkk~ 2 P k . (31) 

Jo 

The correlation parameter is 

B=(t 1 t 2 )/(t 2 ) = dkk- 2S -j^P k / J dkk~ 2 P kl (32) 

where L = 2H ~ 1 [1 - (1 + z)' 1 / 2 ] is the comoving distance to the last scattering surface 
at redshift z. (The small angle approximation $ <C 1 is used throughout.) From (18) the 
coherence angle is then 



"OO 

e 2 t =3L~ 2 j dkk- 2 P k j j dkP kl (33) 
or t « l/(kL). 

For the ray deflection variables, the physical basis is gravitational lensing due to the 
density fluctuations along the line of sight. Following Linder (1990b) the geodesic equation 
gives 

pL px 

9 = 2L 



1 [ dx [ dx'V±(j>, (34) 
Jo Jo 



i.e. the deflection is caused by gradients in the gravitational potential. The equation of 
geodesic deviation governs the relative deflection, or deviation, of neighboring rays: 



3 = 02-0! = 2L" 1 f dx [ dx' x ($ • V±) V±0. 
Jo Jo 



(35) 



Converting the gravitational potential to the density power spectrum as before, by 
^ = (2tt)- 3 / 2 fd 3 k(l) k e i ^ = (3/2)(2tv)- 3 / 2 H 2 f d 3 kk- 2 d k e^ x , one then carries out the 
derivatives and integrations in (34) and (35). To calculate a variance one does the angle 
and ensemble averages as represented by the angle brackets to obtain for the variance of a 
single ray deflection 

poo 

(36) 



poo 

a 2 = (9 2 ) = (3/2n)H i L / dkk~ x P k , 

Jo 
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and for the variance of the ray bundle deviation 



S 2 = (D 2 ) = 2s 2 = 2a$ (1 - b e ), (37) 
in agreement with (12). The deflection correlation parameter comes from (35)-(37): 



b e al = J Q dk k ~ lp k J d 2 n f(ti, O', kL) 

poo 

= (3/2ir)H*L / dk k' 1 [1 - (kL$) 2 /40 + <D(kL$) 4 ] P k , 
Jo 



(38) 



where the angular variables fi, O' of k relative to the two lines of sight are related by a 
rotation of the axes by angle $. For example the direction cosine of the wave vector with 
respect to the second line of sight is \l' = fi cos$ + \J\ — /J, 2 sin $ cos (3 with \i the cosine 
with respect to the first line of sight and (3 the azimuthal angle in the first frame. 

Unfortunately the general evaluation of (38) is not analytically amenable but the 
second line of the equation shows the behavior for angles $ <C (kL) -1 . This is sufficient 
to find the coherence angle and other interesting physical parameters like the gravitational 
lensing ray crossing variable V 2 = limq,^ Q (D 2 ) /$ 2 = (3/407r)i7o-^ 3 / dkkP k . When V 2 > 
1 then the Jacobian of the sky-image map goes singular, rays cross, and the image appears 
grossly distorted or multiple. The coherence angle of the lensing deflections is defined 
analogously to (33), 

poo poo 

Q 2 e = 20L- 2 J dkk- x P k j J dkkP k , (39) 

or 9 w l/(kL). 

For the cross correlation parameter b one takes t = <f>/3 and (35) for D, both involv- 
ing Fourier transforms of 4>k and hence 5k, and carries out the same ensemble average 
calculational process. This gives 



poo 

bsa = (D\\t) = -(3/4Tv 2 )H^a / dk k~ 2 P k [1 - (l/2)k 2 L 2 ^ 2 + 0{kL^)% 

Jo 



(40) 



with a = 3> • ip and again the general behavior unamenable. In the limit of small field 
separation s ~ $ (see equations 37, 38) so 



b^-iUO/n^aL- 3 / 2 



poo poo 

J dkk- 2 P k j J dkkP k 



1/2 



(41) 



For the larger angle behavior one can turn to the more physical intuitive (though not more 
calculationally tractable) expression 

b 2 = [(D 2 t 2 ) - (D 2 )(t 2 )] /(2s 2 a 2 ) = {D 2 t 2 )/(D 2 ){t 2 ) - 1, (42) 

which follows from the joint probability function. This shows that indeed b vanishes as the 
deflections and temperature fluctuations become uncorrelated. 
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To obtain specific quantitative results for the observed microwave background tem- 
perature anisotropy correlation function requires becoming more model dependent, choos- 
ing a particular density power spectrum. We illustrate the orders of magnitude for the 
variables and the behavior of the resulting correlation function by adopting a cold dark 
matter spectrum normalized to the detected microwave background quadrupole anisotropy, 
P k = 1.8 x 10- 9 (a 2 /2.2 x lO" 5 ) 2 ^" 4 k/[l + hk + l 2 k 3 / 2 + Z 3 /c 2 ] 2 , where h, Z 2 , h are length 
constants (Davis et al. 1985). The results are a = 1.1 x 10 -5 , <dt = 14', ue = 3', 0# = 5', 
and b = —2 x 10 -4 . This value for b is at $ = but its maximum at finite $ is only a few 
times larger. In terms of approximating the coherence angles by (ZcL) -1 , the characteristic 
wavenumber is k (15 h~ 1 Mpc)~ 1 . 

5. Conclusion 

We have investigated the behavior of the observed CMB temperature anisotropy correla- 
tion function in a universe with density inhomogeneities causing a gaussian light deflection 
field, both correlated with and independent of the perturbations giving rise to the tem- 
perature fluctuations. The ratio of the correlated effect to the independent one is of order 
the cross correlation parameter squared, b 2 . While one might expect that the deflection 
and temperature be strongly coupled, 6^1 (cf. equation 42), this is found not to be true 
here. In the case of gravitational lensing of photons from the last scattering surface the 
correlation influence on CMB anisotropies is negligible because one loses the path length 
"resonance" of the deflections and temperature fluctuations together tracing the gravita- 
tional potential field along the entire line of sight. Such a dependence would increase the 
correlation parameter by a factor of « kL ps 10 2 5 — 10 3 5 , making the effect significant. 
Due to the quasistatic nature of the potential along the path, however, this does not give 
rise to observable temperature fluctuations. Exceptions occur only in localized events such 
as cluster formation (Rees & Sciama 1968) or hot cluster cores (Sunyaev & Zel'dovich 
1970), which do not offer the opportunity for path resonance or gaussian statistics. 

Although the correlation of the two density perturbation derived variables considered 
here can be neglected, the formalism of calculating probability distributions and covari- 
ances of gaussian fluctuations is powerful. Insights garnered here into the methods and 
relations may prove useful in other cosmological applications such as analyzing large scale 
velocity flows. 

Acknowledgments. I thank Jacob Sharpe for a careful reading of the paper and the referee, 
Jim Bartlett, for improvements in clarity. 
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Figure caption 

Fig. 1. The geometry of the intrinsic field separation vector if), the deflection vectors 9 
for the two lines of sight, and the observed separation vector $ is shown in the plane of 
the sky. 
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